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1. Introduction 



The study of the chemical evolution of gas and dust from pre-stellar dense 
^j^l cores to circumstellar disks around young stars forms an essential part of 

O ' understanding star- and planet formation. Throughout the collapse- and 

protostellar phases, simple and complex molecules are formed, many of 
^ I which deplete onto cold grains and are eventually incorporated into the icy 

planetesimals of new solar systems (see Figure 1). Tracing this chemical 
evolution provides a wealth of information, not only about the chemical 
rN I processing in primitive solar nebulae, but also about physical processes 

j^ ■ which occur in the immediate surroundings of young stellar objects (YSOs). 

Interstellar chemistry has been an active field of research for more than 
60 years, ever since the detection of the first molecules in 1937-1941. Nearly 
120 different gas-phase species have been identified (see Table 1), not in- 
cluding isotopic varieties, with abundances down to 10~^^ with respect to 
H2. The majority of the species have been detected through their rotational 
transitions at millimeter wavelengths. Most of the early observations were 
performed with typical beam sizes of 1', corresponding to linear scales of 
nearly 10,000 AU (0.04 pc) in the nearest star-forming regions in Taurus and 
Ophiuchus. With the advent of large single-dish submillimeter telescopes 
and millimeter interferometers, it has become possible to study the dense 
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Figure 1, A schematic view of the characteristic molecules at different stages of low-mass 
star formation, (a) Dark cloud cores, where radicals and unsaturated carbon chains are 
prominent; (b) Collapse stage, where high levels of depletion of molecules have been in- 
ferred; (c) Deeply embedded YSO phase, where heating and supersonic outflows result in 
evaporation of ices and create a variety of chemical environments; (d) Young T Tauri star 
with a residual protoplanetary accretion disk containing gas-phase molecules and ices; (e) 
Mature planetary system with icy bodies such as comets (Figure by M.R. Hogerheijde, 
after Shu et al. 1987; from: van Dishoeck & Blake 1998). 



envelopes around YSOs more directly on scales of ~ 2"-20" (300-3000 AU 
at 150 pc) (Blake 1997) (see Figure 2). 

Although much effort has focused on gas-phase chemistry, water ice was 
identified already in 1973 through its vibrational absorption toward bright 
infrared sources. The enormous advances in infrared instrumentation in 
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the last decade, both from the ground and from space with the Infrared 
Space Observatory (ISO), have resulted in the detection of several other 
solid-state species and allow a complete inventory to be made (Tielens & 
Whittet 1997). The development of realistic models which include both 
gas-phase and grain-surface chemistry provides a challenge to theorists. 

The combination of new observations and models has lead to the follow- 
ing scenario for the chemistry during star formation (see Figure 1). In the 
cold molecular cores prior to star formation, the chemistry is dominated 
by low-temperature gas-phase reactions leading to the enhanced formation 
of small radicals and unsaturated molecules. Long carbon-chains are pro- 
duced if the gas is initially atomic-carbon rich. During the cold pre-stellar 
and collapse phases, many molecules accrete onto the grains and form an 
icy mantle. Here surface chemistry and processing by ultraviolet photons. 
X-rays and cosmic rays can modify the composition. Once the new star 
starts to warm the surrounding envelope, the ices are heated and molecules 
evaporate back into the gas phase, probably in a sequence according to 
their sublimation temperatures. In addition, the outflows from the young 
star penetrate the envelope, creating high temperature shocks and lower 
temperature turbulent regions in which both volatile and refractory mate- 
rial containing silicon can be returned. These freshly evaporated molecules 
then drive a rich chemistry in the 'hot cores' for a period of 10^ yr. Finally, 
the envelope is dispersed by winds and/or ultraviolet photons, leading to 
the appearance of photon-dominated regions. 

Part of the gas-phase and icy material is incorporated into the circum- 
stellar disks, where they survive for a few xlO^ yr before being dispersed 
or assembled into new planetary bodies. Indeed, the connection between 
interstellar chemistry and solar system material has been greatly strength- 
ened by observations of comets Hyakutake and Hale-Bopp, which reveal a 
remarkable similarity between the composition of interstellar and cometary 
ices (e.g., Ehrenfreund et al. 1997, Irvine et al. 1999). 

Molecular spectroscopy is also a very powerful tool to probe the enor- 
mous range in temperatures from 10 to 10,000 K and densities from 10*^ to 
10^^ cm~^ that accompany the star-formation process. In addition, molec- 
ular lines provide the only means to trace the velocity fields in the highly- 
extincted circumstellar environment. In fact, the two analyses go hand in 
hand. On the one hand, an accurate physical model is a prerequisite for 
deriving reliable molecular abundances. On the other hand, it is important 
to know which molecule traces which physical component if its emission is 
to be used as a probe of the physical parameters. 

This chapter first discusses the basic physical and chemical processes 
that play a role during star formation. Subsequently, the observational 
methods and problems in deriving of abundances are discussed. Finally, 
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Table 1 Identified interstellar and circumstellar molecules 



Simple Hyd 


rides, Oxides, 


Sulfides, Halogens and related molecules 


H2 (IR) 


CO 


NH3 


CS 


NaCl* 


HCl 


SiO 


SiH* (IR) 


SiS 


Aicr 


H2O 


S02 


C2 (IR) 


H2S 


KCl* 


N2O 


ocs 


CH4 (IR) 


FN 


AlF* 


HF 










Nitriles and Acetylene derivatives 






C3 (IR.UV) 


HCN 


CH3CN 


HNC 


C2HJ (IR) 


C5 (IR) 
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CH3C3N 
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CH3C5N ? 
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C4Si-' 
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HCiiN 


CH3CH2CN 
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HC2CHO 


CH2CHCN 






Aldehydes, 


Alcohols, Eth< 


Brs, Ketones, Amides and relat* 


Bd molecules 


H2C0 


CH3OH 
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CH2NH 


CH2CC 
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CH3COOH 




Cyclic Molecules 
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C-C3H 
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ons 
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H2COH+ 


S0+ 






Radicals 
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From Ohishi (1997) with detections as of October 1998 added; Species denoted with * have 
only been detected in the circumstellar envelope of carbon-rich stars. Most molecules have 
been detected at radio and millimeter wavelengths, unless otherwise indicated (IR, VIS or 
UV). 



recent results on observations and models of the chemistry in the vari- 
ous star-formation phases are presented. This discussion focuses on those 
species whose abundances are particularly enhanced or decreased at a cer- 
tain phase. These species often have only minor abundances in terms of 
overall composition — less than 10"^ with respect of hydrogen — , but pro- 
vide important 'signposts' or 'clocks' of the evolutionary state of an object. 
This chapter is based on the recent reviews by van Dishoeck (1998a) on 
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Figure 2. JCMT submilllimeter spectra of the deeply embedded low-mass YSO IRAS 
16293—2422, showing lines of characteristic molecules such as H2CO, CH3OH and deuter- 
ated species (from: van Dishoeck et al. 1995). 



the chemistry in diffuse and dense clouds, and by van Dishoeck & Blake 
(1998), van Dishoeck (1998b) and Langer et al. (1999) on the chemistry 
in star-forming regions. Other reviews on this subject include Hartquist et 
al. (1998), Irvine (1998) and van Dishoeck et al. (1993). For an excellent 
overview of the physics of molecular clouds, see Genzel (1992). 



2. Gas-phase chemistry 

2.1. BASIC MOLECULAR PROCESSES 

The networks developed to describe the gas-phase chemistry in star-forming 
regions contain thousands of reactions between several hundreds of species 
(Millar et al. 1997a, Lee et al. 1996a). These thousands of reactions involve 
only a handful of different basic molecular processes, however, which have 
been described in detail in reviews by Dalgarno (1987) and van Dishoeck 
(1988). At the low densities in the interstellar gas, three-body processes 
are unimportant, so that only two-body reactions need to be considered. 
The rate of a reaction between species X and Y is given by A;n(X)n(Y) in 
cm~^ s"^, where k is the reaction rate coefficient in cm'^ s^"*^ and n is the 
concentration in cm~^. 

There are two basic mechanisms by which molecular bonds can be 
formed (see Table 2). The first process is radiative association of gas-phase 
atoms or molecules, in which the new molecule is stabilized by the emis- 
sion of a photon. The second mechanism is the formation of molecules on 
the surfaces of grains, in which the grain carries off the released energy 
corresponding to the molecular bond. The rate coefficients for the forma- 
tion processes range from less than 10^^^ up to 10^^^ cm^ s^^ and have 
considerable uncertainties, often up to an order of magnitude. 

Molecules are readily destroyed by the absorption of ultraviolet pho- 



EWINE F. VAN DISHOECK AND MICHIEL R. HOGERHEIJDE 



TABLE 2. Some basic microscopic processes 



Process 



Name 



Chemical Processes 



X + Y ^ XY 
X + Y:grain — 



- Yiv 
XY 



gram 



XY + hi/ - 

XY+ + e - 
X+ + YZ 
X+ + YZ 
X + YZ - 



-> X + Y 
-+X + Y 
-^ XY+ + Z 
-^ X + YZ+ 
XY + Z 



Heating and Cooling Processes 

grain^ or PAH"'" + e* 



H++e* 



grain or PAH + \iv 

H2 + cosmic ray -^ 

CO(J) + coll. ^ CO(J*) -^ CO(J') + hi/ 

OC^Pa) + coU. ^ 0(''Pi) ^ 0(^P2) + h// 



c+rPi/2) + C0IL - 

gas + grain -^ gas' 



C+(^P3/2) 

+ grain' 



1/2, 



+ hi/ 



Radiative association 
Grain-surface formation 
Photodissociation 
Dissociative recombination 
Ion-molecule reaction 
Charge-transfer reaction 
Neutral-neutral reaction 



Photoelectric heating 
Cosmic ray heating 
CO line cooling 
[O I] line cooling 
[C H] line cooling 
Gas-grain heating or cooling 



tons. This process is very effective at the edges of dark clouds and in clouds 
located close to young stars (so-called photon-dominated regions (PDRs), 
Hollenbach & Tielens 1997, 1999). For molecules such as H2 and CO, the 
photodissociation can take place only at very short wavelengths between 
912 and 1100 A, whereas most other species are dissociated by radiation out 
to 3000 A. The ultraviolet photons can also ionize atoms with ionization 
potentials less than 13.6 eV (e.g., C, S, Si, Fe, Mg, ...), thereby increasing 
the electron abundance inside the cloud. In the unshielded interstellar radi- 
ation field such as given by Draine (1978), typical rates are 10~^-10"^'^ s^^, 
corresponding to lifetimes against photodissociation or photoionization of 
only lO^-lO^ yr. 

Inside a cloud, the ultraviolet radiation is reduced because of absorption 
and scattering by grains (e.g., Roberge et al. 1991). Deep inside dark clouds 
(> 5 mag), little of the ambient radiation penetrates, but a weak ultravi- 
olet field is maintained by cosmic-ray induced photons, resulting from the 
excitation of H2 by secondary electrons produced by cosmic ray ionization 
of H2 (Gredel et al. 1989). For a typical cosmic ray ionization rate C of a few 
xlO~^^ s~^, the resulting photorates are four to five orders of magnitude 
lower than those in the unshielded interstellar radiation field at the edge. 
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Figure 3. Initial steps in the gas-phase carbon chemistry in diffuse and dark clouds. 



Molecular ions are efficiently destroyed by the process of dissociative 
recombination. Rate coefficients in 10-30 K gas are typically 10~^ to 10~^ 
cm^ s~^, with the lower values representative for Hj]" — a key species in 
the chemistry (Dalgarno 1994, Larsson 1997). A major uncertainty in the 
models is the branching ratio to the various products. Complete experi- 
mental data are only just becoming available, and recent results indicate 
that three-body product channels (e.g., H3O+ + e — > OH + H -|- H or O 
+ H2 + H) have a much larger probability than thought previously (e.g., 
Vejby-Christensen et al. 1997). 

In dense clouds, destruction of neutral molecules can also occur through 
chemical reactions, in particular with He^ (e-g-, He"*" + CO — > He + C^ 
+ O). Collisional dissociation is only important in regions of very high 
temperature (> 3000 K) and density, such as found in shocks. 

Once molecular bonds have been formed, they can be rearranged by 
chemical reactions leading to more complex species. The first generation of 
models (e.g., Herbst &: Klemperer 1973, Prasad & Huntress 1980) included 
primarily ion-molecule reactions. The majority of these reactions are known 
to be fast at low temperatures owing to the long-range attraction between 
the ion and the molecule. It has recently been realized that neutral-neutral 
reactions involving radicals and/or atoms can also have rate coefficients 
as large as ~ lO^^^'-lO"^ cm^ s~^ at low temperatures (e.g.. Smith 1997, 
Kaiser et al. 1997). 
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2.2. GAS-PHASE NETWORKS 

The gas-phase chemical networks describing the most important forma- 
tion routes of various molecules based on the above processes have been 
discussed extensively in the literature, see reviews by, e.g., Millar (1990), 
Herbst (1995) and van Dishoeck (1998a). The abundances of the elements 
play a key role in the chemistry, especially the absolute and relative amounts 
of gas-phase carbon and oxygen in atomic or molecular form. It has recently 
become clear that the interstellar abundances of most elements are ~ 30% 
lower than the solar abundances (see Meyer 1997 for an overview). The 
fraction of heavy elements in the gas phase is large (at least 30%) in diffuse 
clouds, except for the most refractory elements (Si, Fe, ..), which are 99% 
in solid form. For dense clouds, the depletions of species like C, N, O and S 
are expected to be more significant, but are ill constrained observationally. 

Because hydrogen is so much more abundant than any other element, 
reactions with H and H2 dominate the networks if they are exothermic. 
This is only the case for small ions at low temperatures; most reactions of 
neutrals and large ions with H or H2 have substantial barriers. H2 itself is 
produced predominantly on the surfaces of grains (Hollenbach &: Salpeter 
1971), but once formed it is essential for the gas-phase chemistry of other 
species. Ionization of H2 by cosmic rays at a rate C leads to Hg , which 
reacts with H2 to form the stable H3 ion. This ion plays a pivotal role in 
the subsequent ion-molecule chemistry through proton transfer. Its recent 
detection in interstellar space provides strong observational support for the 
gas-phase chemical networks (Geballe & Oka 1996, McCall et al. 1998). 

The production of complex hydrocarbons in dense clouds occurs via 
three types of pathways (Herbst 1995): (i) carbon insertion reactions (e.g. 
C+ + CH4 -^ C2H+ + H2 or C + C2H2 -^ C3H + H); (ii) condensation 
reactions (e.g. CH+ + CH4 -^ C2H+ + H2 or C2H + C2H2 -^ C4H2 + H); 
and (iii) radiative association reactions (e.g. C"*" + C„ ^ ^n+i + ^^) (^^^ 
also Figure 3). In general, carbon insertion with C"*" is thought to be the 
dominant route. Since this leads to loss of one hydrogen and because the 
larger ions C„II^ do not react rapidly with H2, low-temperature gas-phase 
chemistry produces strongly unsaturated hydrocarbons, in agreement with 
observations of dark clouds. The necessary C''' ions are produced in small 
amounts through destruction of CO by He"*". Reactions with C may be com- 
petitive if they are as rapid as suggested by recent laboratory experiments. 
Once carbon is locked up in the very stable CO molecule, the formation of 
more complex hydrocarbons ceases. 

At high temperatures of ~ 200-2000 K such as encountered in the 'hot 
core' regions surrounding massive young stars and in shocks, gas-phase 
reactions with H and II2 become significant. Specifically, the reactions O + 
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H2 -^ OH + H - 4480 K and OH + H2 -^ H2O + H - 2100 K 
start to proceed at temperatures of a few hundred K (Ceccarelli et al. 
1996, Charnley 1997), and drive most of the available gas-phase oxygen 
not locked up in CO into H2O. Similar reactions of S with H2 convert the 
gas-phase sulfur into H2S. However, the back reactions of OH and H2O 
with H have comparable energy barriers, so that the balance between O, 
OH and H2O depends also on the H/H2 ratio in the warm gas. 

When most of the oxygen has been driven into H2O at high tempera- 
tures, little O, OH and O2 are available in the gas, preventing the formation 
of, for example, S02- Since O and O2 also destroy reactive species such as 
CN, the CN abundance is enhanced at higher temperatures leading to en- 
hanced formation of HC3N through reaction with C2H2. 

2.3. GAS-PHASE MODELS 

2.3.1. Depth- dependent vs. time- dependent models 

The calculation of the abundances in star-forming regions requires a phys- 
ical model in which the temperature, density, radiation field etc. are speci- 
fied as functions of position and/or time. In general, two different classes of 
models are considered: (i) Steady-state, depth-dependent models, in which 
the physical parameters and molecular abundances do not change with 
time, but are functions of depth into the cloud. Models of the translucent 
outer envelopes of clouds (e.g., van Dishoeck 1998a) and of dense ultravi- 
olet photon- or X-ray dominated regions (PDRs or XDRs, see e.g., Hol- 
lenbach & Tielens 1997, Sternberg et al. 1997, Maloney et al. 1996) fall 
in this category, (ii) Time-dependent, depth-independent models, in which 
the concentrations are computed as functions of time at a single position 
deep into the cloud. Models of dark pre-stellar cores (e.g., Lee et al. 1996a, 
Millar et al. 1997b), collapsing envelopes (e.g., Bergin & Langer 1997) and 
hot cores near massive young stars (e.g., Charnley et al. 1992) fall into this 
category. The time scale for reaching chemical equilibrium depends on the 
density, temperature and ionization fraction and is typically 10^-10^ yr. 

In both cases, the parameters that enter the models are (1) the elemental 
abundances of C, O, N, S, metals...; and (2) the cosmic ray ionization rate C 
in s~^. In the steady-state, depth-dependent models, additional parameters 
are (3) the geometry (e.g. plane-parallel, spherical, ...); (4) the density 
'^H=^(H) -|- 2n(H2) as a function of position; (5) the incident radiation 
field, specified by a factor Ijjv times the standard interstellar radiation 
field as given by, e.g., Draine (1978); and (6) the grain parameters, i.e., 
the extinction curve, albedo and scattering function. The temperatures of 
the gas and dust as functions of position in the cloud can be obtained self- 
consistently from the balance of the heating and cooling processes listed 
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in Table 2. Alternatively, they can be constrained from observations and 
provided as additional input parameters (see §4.2). 

In the time-dependent models, the additional parameters besides (1) 
and (2) are: (3) the density, usually taken to be constant with time (so- 
called pseudo time-dependent models); (4) the visual extinction Ay at the 
position in the cloud, usually taken to be so large that external photons can 
be neglected; and (5) the initial abundances of the various species at t=0, 
often taken to be atomic except for H2. The temperature can be obtained 
from the thermal balance, but is almost always set at 10 K for both the gas 
and dust, typical of a dark cloud shielded from ultraviolet radiation and 
heated by cosmic rays only. In comparison with observations, the ratios of 
the local concentrations (in cm~^) are taken to be equal to the ratios of 
the column densities integrated over depth (in cm^^). This assumption is 
accurate for molecules whose abundances peak in the center of the cloud, 
but not necessarily for species such as radicals whose abundances peak in 
the outer part of the envelope (see Figure 5 of van Dishoeck 1998a). 

2.3.2. The C+ ^ C -^ CO transition 

The principal chemical characteristics of the depth- and time-dependent 
models are governed by the transition of carbon from atomic to molecular 
form. In the depth-dependent case, C"^ recombines to C around Ay ~ 1 
mag, followed by the conversion to CO around Ay ~ 2 mag. The CO 
photodissociation rate as a function of depth is crucial in this transition 
(e.g., van Dishoeck &: Black 1988, Lee et al. 1996b). The chemistry of other 
species follows the C"^ — > CO transition. At the edge, only the simplest 
diatomic molecules are found. Around Ay = 1-2 mag, the presence of both 
C and C+ and simple hydrides results in an increase in the abundance 
of hydrocarbon molecules such as CN and C2H. Deeper inside the cloud, 
atomic carbon is no longer available, and destruction by O becomes more 
important than photodissociation. Stable species such as CH4, C2H2 and 
HCN become dominant. 

Many of the same features are observed in the pseudo time-dependent 
models, if depth is replaced by time. These models usually start with the 
assumption that dark clouds originate from diffuse gas, so that all species 
except H2 are initially in atomic form, with carbon present as C"*". On a 
time scale of ^ 10^ yr, C"^ recombines to C, which subsequently transforms 
to CO after --^ 10^ yr. Since the presence of atomic carbon is essential to 
build up more complex organic species such as IIC3N, they are abundant 
only at early times, but not at steady-state (see Figure 4). 

The main time-dependent aspects of the other elemental pools — O, N 
and S — , depend on how their chemistry is linked to that of carbon. For 
example, the abundance of CS has a different time behavior than that of 
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Figure 4- Time dependence of the abundances of C^, C and CO (left) and HC3N 
(right) for a cloud with wh = 2 x 10'' cm~^. The latter behavior is typical of all complex 
carbon-bearing molecules (from: Herbst & Leung 1986). 



SO, because CS is formed by reactions with carbon- hydrides, whereas SO 
is destroyed by atomic carbon. This leads to two classes of molecules: (i) 
species like CN, HCN, CS and complex carbon chains, whose production 
is linked to the carbon network and which have larger abundances at early 
times; and (ii) molecules such as N2, NH3, N2H"'" and SO that are inde- 
pendent or destructively linked to the carbon chemistry and which exhibit 
higher abundances at equilibrium. 



2.4. FRACTIONATION 



At the low temperatures of pre-stellar cores and the outer envelopes of 
YSOs, significant enhancement of minor isotopes can occur as a result of so- 
called 'fractionation'. The most relevant case for star- forming regions is the 
deuterium fractionation, initiated by the isotope exchange reactions H;^ -|- 
HD -^ H2D+ + H2 and CH+ -h HD ^ CH2D+ + H2 (Wootten 1987, Millar 
et al. 1989). These reactions occur preferentially in the forward direction at 
low temperatures, because the heavier H2D+ and CH2D^ are slightly more 
stable than Hg" and CHg" due to their lower zero-point vibrational energy. 
H2D^ and CH2D''' subsequently transfer the deuterium to, for example, 
CO or N to form DCO+ or DCNH+. The latter ion leads to DCN through 
dissociative recombination. The observed DCO"'"/HCO"'" and DCN/HCN 
abundance ratios are factors of 100 to >1000 larger than the overall [D]/[H] 
ratio of ~ 1.6 X lO^^. 
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2.5. SUCCESSES AND PROBLEMS 

Pure gas-phase chemistry models have been remarkably successful in ex- 
plaining many observed features of interstellar chemistry (Herbst 1995, 
Turner et al. 1998). These include the abundances of simple species in dif- 
fuse and translucent clouds; the presence of Hg and the related high abun- 
dances of protonated species such as HCO^ and N2H^ in dark clouds; the 
high abundances of unsaturated and metastable molecules in dark clouds 
such as /— and c— C3H or HNC; and the large isotopic fractionation. 

Most of the problems with gas-phase models stem from uncertainties in 
the basic reaction rates. For example, if reactions of bare carbon chains C^ 
with O or N are rapid, the observed abundances of carbon-chain molecules 
like HC7N cannot be reproduced, even at early times. Another complica- 
tion is that time-dependent models have more than one solution in certain 
regions of parameter space (e.g., le Bourlot et al. 1995, Lee et al. 1998). It is 
not yet clear, however, what the astrophysical consequences of this 'bista- 
bility' phenomenon are. Finally, gas-phase models of dark clouds predict 
that a substantial fraction of the oxygen is driven into O2 at low temper- 
atures. Recent limits on the O2/CO ratio of less than 0.03 are difficult to 
accomodate in such models (Olofsson et al. 1998, Marechal et al. 1997, 
Melnick et al. 1998). 

3. Grain-surface chemistry 

3.1. BASIC SURFACE PROCESSES 

The chemistry on the surfaces of interstellar grains has received ample dis- 
cussion in the literature (e.g., Tielens & Allamandola 1987, Herbst 1993, 
Tielens &: Charnley 1997, Langer et al. 1999). Four different steps can be 
distinguished: (i) accretion; (ii) diffusion; (iii) reaction; and (iv) ejection or 
evaporation. The time scale for a molecule to collide with a grain and ac- 
crete is given by tacc ~ 2 x lO^/nnys yr where ys is the sticking probability 
which is thought to lie between 0.1 and 1.0 for most species (see Williams 
1993 for a review). Thus, in envelopes around YSOs with densities of order 
10^ cm~^ or larger, the time scale for depletion is shorter than the col- 
lapse time (Walmsley 1991). The timescales for diffusion and evaporation 
increase exponentially with the binding energies of the species to the grain, 
which in turn scale with its mass in the case of physical adsorption. Under 
most interstellar circumstances, the rate-limiting step is the accretion of 
new species rather than the diffusion over the surface. 

Because atomic hydrogen is so abundant and mobile, grain-surface chem- 
istry leads primarily to hydrogenated species such as H2O, NH3 and CH4. 
Although many of the surface reactions have energy barriers, they can still 
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Figure 5. Grain-surface chemistry routes involving CO. Solid rectangular boxes contain 
molecules which have been identified in interstellar ices, whereas dashed boxes indicate 
molecules detected in the gas phase in hot cores (figure by Tielens, based on Tielens & 
Hagen f982). 



proceed because of the relatively long time scale (~ 1 day) available for 
reaction before another species is delivered from the gas to the grain. 

The dominant molecule accreted from the gas is CO, and it there- 
fore forms a key in the surface chemistry. Hydrogenation of CO leads 
to HCO, H2CO and eventually CH3OH (see Figure 5). The intermediate 
HCO (formyl) and CH3OH (methoxy) radicals can also react with vari- 
ous other accreted species, leading to more complex organic molecules such 
as CH3CHO, HCOOH and HCONH2. The identification of many of these 
species in molecular clouds, directly in ices or indirectly in hot cores where 
the ices have evaporated off the grains, provides support for this scheme. 

Hydrogenation of molecules is particularly effective in the outer regions 
of YSO envelopes where the atomic hydrogen abundance is high. At den- 
sities greater than ~ 10^ cm~^ shielded from ultraviolet radiation, atomic 
oxygen becomes more abundant than atomic hydrogen, leading to oxida- 
tion of various accreted molecules. For example, CO can react with O to 
form CO2, and O2 with O leads to O3. 

Deuterium fractionation of molecules can also occur through grain- 
surface chemistry, because the gas-phase atomic D concentration inside 
dense clouds is enhanced compared with atomic H through dissociative re- 
combination of DCO"*" (Tielens 1983). The atomic D can then react with 
molecules on the surface to form deuterated species. Doubly-deuterated 
molecules such as D2CO and NHD2, which have been detected in the en- 
velopes around low- and high-mass stars, are thought to be formed by this 
route (e.g., Ceccarelli et al. 1998). 
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TABLE 3. Interstellar ice composi- 
tion and evaporation temperatures 



Species 


Abundance" 


Tev' (K) 


H2O 


100 


90 


CO 


16 


16 


CO2 


24 


50 


CH4 


2 


18 


CH3OH 


4 


80 


XCN 


2 


7 


OCS 




85 


NH3 


10= 


55 


N2 




13 


O2 


OC* 


13 



" Abundances are relative to H2O=100, and refer to NGC 7538 IRS9 (Whittet et al. 

1996). Abundances are variable from source to source due to differential out gassing and 

thermal processing. The H2O abundance is typically 10~ with respect to H2. 

'' Evaporation temperatures for pure ices under interstellar conditions; see text for mixed 

ices. 

'^ Tentative identification by Lacy et al. (1998). 

'^ Vandenbussche et al. (1999). 



3.2. THERMAL PROCESSING AND EVAPORATION 

In star-forming regions, the ice mantles formed in the cold pre-stellar and 
collapse phases can be heated by the young star. This leads to restructuring 
of the ice matrix and outgassing when the temperature is near its subli- 
mation point. Table 3 summarizes the sublimation temperatures of various 
species under interstellar conditions (Sandford &: Allamandola 1990, 1993); 
these are typically lower than found in laboratory experiments because of 
the lower pressures in space. Interstellar ices are not pure, but in mixed 
ices, the evaporation of each component is largely determined by its own 
sublimation behavior, unless the abundance of one of the species is less 
than rsj 5%. Thus, an ice mixture of H20/CO/C02=2/l/l shows release 
of CO around 20 K, CO2 around 50 K and H2O around 90 K. However, 
a mixture of H2O/CH4=100/l shows evaporation of both H2O and CH4 
around 90 K. Only a small fraction of the CO or CO2 stays behind, and is 
released when the whole H2O ice rearranges from amorphous to crystalline 
form. 

Ices containing CH3OH show particularly interesting behavior upon 
heating, because CH3OH has a sublimation temperature above that of the 
phase transformation of H2O. Thus, H2O/CH3OH ices heated to ~ 80- 
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90 K will segregate into rather pure H2O- and CH3OH domains (Blake et 
al. 1991). Recent experiments indicate that CO2/CH3OH mixtures show 
similar segregation behavior (Ehrenfreund et al. 1998), and evidence for 
this process is seen in the observed profiles of interstellar solid CO2 toward 
massive YSOs (Boogert et al. 1999, Gerakines et al. 1999). 



3.3. NON-THERMAL DESORPTION 

Thermal evaporation is only efficient at dust temperatures T^nst higher than 
20 K. Some desorption can also occur by non-thermal processes in the cold 
pre-stellar cores and outer envelopes with Tjust ~ 10 K. These processes 
have been reviewed most recently by Schutte (1996), and include cosmic 
ray spot heating, heating due to the energy liberated by the formation 
of molecules on grains, heating due to grain-grain collisions, and explo- 
sive heating due to exothermic reactions between stored radicals on grains. 
The desorption time scales depend strongly on the binding energies of the 
molecules at the surfaces and the surfaces involved (silicates, carbonaceous 
material, ices). All of them are particularly effective for non-polar ices con- 
taining CO, O2 and N2, and less for H20-rich ices which contain strong 
hydrogen bonds. 



3.4. ENERGETIC PROCESSING 

The cosmic-ray induced ultraviolet photons and the radiation produced by 
the young star can further process the ices. Ultraviolet photolysis of H2O, 
CO, NH3 and CH4 produces radicals which can react with each other and 
with the parent molecules to form more complex molecules like H2CO, 
HCOOH and HCOCH3 (e.g., Bernstein et al. 1995, Gerakines et al. 1996). 
Much of the outcome of this chemistry is very similar to that of the hydro- 
genation of CO, and further quantitative studies are needed to assess their 
relative importance. 

Photolysis of ice mixtures containing H2O, CO, NH3 and CH3OH of- 
ten leaves a non-volatile organic residue containing a variety of complex 
molecules (e.g., Schutte et al. 1993, Bernstein et al. 1995). Some of these 
chemical effects are also found in experiments in which ices are bombarded 
with highly-energetic charged particles, analogous to cosmic rays or X-rays 
(e.g., Moore et al. 1983, Kaiser k Roessler 1997, Teixeira et al. 1998). The 
photochemical and high-energy bombardment processes are often referred 
to as 'energetic' processing in the literature, without discrimination. 
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3.5. POLAR AND APOLAR ICES 

The ices in interstellar clouds are generally not homogeneous, but consist 
of different layers or domains (e.g., Tielens et al. 1991, Ehrenfreund et 
al. 1998). Two broad classes of ices can be distinguished. Polar ices are 
dominated by H2O ice and contain minor amounts of CO2, CO, CH4 and 
CH3OH. Apolar ices are dominated by species such as CO and perhaps some 
O2 and N2, but contain very little H2O. The two phases can be distinguished 
by their line profiles, since the force constants of CO molecules embedded 
in an H20-matrix will be slightly different from those of CO surrounded by 
other CO molecules. 

The different ice phases can arise from a combination of at least two pro- 
cesses, caused by the density and temperature gradient in YSO envelopes 
(see Schutte 1999 for review). First, the density gradient results in a steep 
gradient in the gas-phase H/CO ratio. At low densities or at the edge of 
the cloud, the high concentration of H results in polar ices. Deep inside the 
envelope at high densities, mostly apolar ices of accreted CO, O2 and N2 
form. Layered ices can be produced in a collapsing envelope with the polar 
ices condensing first and the non-polar species forming a volatile 'crust'. 
Second, the evaporation processes can shape the composition of the ice 
mantles, because the desorption mechanisms are much more efficient for 
volatile non-polar species like CO than for non-volatile material like H2O 
and CH3OH. This 'distillation' effect decreases the apolar ices compared 
with the polar ices at higher temperatures. 

3.6. GAS-GRAIN MODELS 

Based on the above processes, two different chemical regimes can be distin- 
guished in models which take both gas-phase and grain-surface processes 
into account. In the 'accretion-limited' regime, the diffusion time is much 
shorter than the accretion time so that a species can diffuse on the surface 
until it finds a co-reactant. The chemistry is limited by rate at which new 
species are delivered to the surface. In the 'reaction-limited' regime, the 
opposite holds so that many reactive species are present on a grain surface 
and the reaction is controlled by surface concentrations as well as kinetic 
parameters. Most of the gas-grain chemical models have been formulated in 
the 'reaction-limited' regime using rate equations which mirror those used 
for gas-phase chemistry (e.g., Hasegawa & Herbst 1993). This approach is 
only accurate when a large number of reactive species exist on a single 
grain surface, since only average abundances are calculated. This condition 
is usually not met in interstellar clouds, since the accretion times are long, 
grains are small, and reactions are fast, so that at most one reactive species 
is present on a grain at any time. The surface chemistry is therefore in 
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the accretion-limited regime and can only be properly treated by a Monte- 
Carlo method which determines the likelihood of two such species arriving 
from the gas in succession onto a particular grain in a steady-state model 
(e.g., Tielens & Hagen 1982). Recently, ad-hoc modifications of the rate 
equations have been proposed to correct the shortcomings of the reaction- 
limited approach (Caselli et al. 1998, Shalabiea et al. 1998). 

3.7. SUCCESSES AND PROBLEMS 

The observations of H2 in diffuse clouds and of large amounts of ices in 
dense clouds can only be explained in models in which these molecules are 
formed on the surfaces of interstellar grains and/or accreted from the gas. 
Thus, gas-grain interactions play an essential role in interstellar chemistry. 
The main problem is an accurate, quantitative description. Many basic 
grain-surface reactions with H or O at low temperatures are still poorly un- 
derstood theoretically or ill-constrained by laboratory experiments, partly 
because the details of the surfaces of interstellar grains are still poorly 
characterized. 

4. Determination of molecular abundances 

The comparison of observations and chemical model calculations ultimately 
depends on the translation of spectral- line data into molecular abundances. 
In the following, the observational techniques will be discussed, and the 
methods for constraining the physical structure and subsequently the chem- 
ical abundances will be outlined. 

4.1. OBSERVATIONAL TECHNIQUES 

The energy level structure of a molecule can be decomposed into an elec- 
tronic part Eel, a component due to the vibrational motion of the nuclei 
-Evib; and a component describing the overall rotation of the molecule in 
space Ej-ot- In general, E^i » E'vib >> Ej-ot'- transitions between two elec- 
tronic levels lie typically in the visible/ultraviolet part of the spectrum, 
those between two vibrational levels in the infrared part, and those be- 
tween two rotational levels at (sub-)millimeter wavelengths. Because of the 
high extinction, only infrared and millimeter observations are feasible in 
star-forming regions. 

4.1.1. Rotational line emission at (sub-)miUimeter wavelengths 
Most molecules listed in Table 1 are observed by their rotational transitions 
within the ground electronic and vibrational state (Figure 6). The sensitiv- 
ity of large single-dish (sub-)millimeter telescopes allows the detection of 
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Figure 6. Left: Schematic illustration of infrared absorption line observations of gas and 
dust toward embedded or background sources. The infrared continuum is provided by 
the hot dust at 300-1000 K. Right: Normalized ISO-SWS spectra toward three massive 
young stellar objects embedded in dense molecular clouds. The strong, broad absorption 
at 4.27 /im is due to solid CO2, whereas the characteristic ro- vibrational P- and R-branch 
structure at 4.4-4.9 ^m indicates the presence of warm, gaseous CO along the line of 
sight (van Dishoeck et al. 1998). 



weak lines and rare chemical species and isotopes down to abundances of 
10^^^ with respect to hydrogen. However, not all molecules can be observed 
by this method. Symmetric molecules lack the permanent dipole moment 
required for rotational emission lines, so that important species like H2, 
N2, CH4 and C2H2 cannot be observed in this way. In addition, the Earth's 
atmosphere prevents ground-based observations of low-lying transitions of 
O2, CO2 and H2O. Because of its chemical stability, CO often serves as a 
tracer of H2, assuming either a fiducial CO/H2 abundance ratio or a con- 
version factor between integrated CO intensity and H2 column density (see 
van Dishoeck & Black 1987, van Dishoeck 1998a for reviews). The only 
direct measurement of CO/H2=2.7 x 10"*^ in the warm dense cloud toward 
NGC 2024 IRS2 (Lacy et al. 1994) is a factor of 3 higher than the value of 
8 X 10~^ commonly used in dark clouds (Frerking et al. 1982). SubmiUime- 
ter continuum emission from dust can also act as a standard when values 
are adopted for the dust emissivity, its temperature, and the dust-to-gas 
mass ratio. 

Clouds can be mapped at angular resolutions between 10" (1500 AU 
at 150 pc) and a few arcminutes with single-dish instruments. Even in the 
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nearest star-forming clouds, regions of different physical or chemical char- 
acteristics (envelope, outflow, disk, ...) are likely to fall within one beam, 
complicating the derivation of molecular abundances. Much higher angular 
resolution of the order of one arcsecond is offered by interferometers. How- 
ever, the sensitivity of these instruments is lower than that of single-dish 
telescopes, and their interpretation is complicated because they filter out 
all emission on extended scales defined by the shortest baselines (cf. §8.1). 

A major advantage of heterodyne (sub-) millimeter spectroscopy is the 
very high spectral resolution, typically better than a few tenths of a km s~^ 
and covering the full extent of the lines up to a few hundred km s~^. 
This allows separation of regions with different physical conditions within 
one beam if they have different velocities or line profiles. A drawback of 
rotational emission spectroscopy is that the frequency separation between 
different transitions of the same molecule is usually too large to be obtained 
in a single setting, and multiple observations are needed to constrain the 
excitation, often involving different telescopes or beam sizes. 

4.1.2. Vibrational absorption at infrared wavelengths 

Absorption of infrared emission into (ro-)vibrationally excited levels of 
molecules offers a complementary view of the chemical content (see Fig- 
ure 6). Background stars or embedded objects of sufficient infrared bright- 
ness can usually be found in star-forming regions, often the YSO under 
study itself. Absorption lines probe molecular gas along a pencil beam to 
the infrared source with very high angular resolution, but mapping is not 
possible and different conditions still exist along the line of sight. The ve- 
locity resolution is generally lower, so that the line profiles are often not 
resolved. 

Because of the low spectral resolution and the intrinsic strengths of 
the transitions, the abundance limit of infrared observations is two to four 
orders of magnitude higher than that at millimeter wavelengths. On the 
other hand, symmetric species like C2H2 and CH4 can be observed, be- 
cause the vibrational modes induce a temporary dipole moment. Different 
transitions of a molecule often fall in a single infrared spectrum, providing 
direct constraints on the molecular excitation up to high energy levels. By 
their very nature, absorption lines also directly measure the line opacity, if 
the profiles are resolved. Another major advantage is that infrared studies 
are not limited to the gas phase. Molecules frozen onto dust grains show 
a markedly different line profile, since the rotational degrees of freedom 
are collapsed into a single, broad feature which is slightly shifted from the 
position of that of the gas-phase molecule (Figure 6). This offers the pos- 
sibility to study the relative amounts of a species present in the gas-phase 
and in the condensed phase, and to probe the ice mantles as formation 
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sites of various molecules like CH3OH and H2CO. The Short Wavelength 
Spectrometer (SWS) on ISO is particularly well suited for such studies (see 
§6.2, Figure 6 and 10). 

4.2. CONSTRAINING THE PHYSICAL STRUCTURE 

In order derive molecular column densities or abundances from the obser- 
vations, a good physical model of the region is a prerequisite. Observations 
of different rotational lines provide very useful probes, because they are 
coUisionally excited by H2 for densities of 10^-10 cm^^ and temperatures 
of 5-200 K (Evans 1980, Walmsley 1987, Genzel 1992). 

The statistical equilibrium equations for the level populations rii due to 
radiative and collisional processes are given by 

rii Y^ Aij + ni'Y^{BijJ + CijU^ox) = ^ njAji + ^ nj{BjiJ + CjiUcoi), 

j<i j j>i j 

where the left-hand side represents the processes leading to the loss of 
population of level i and the right-hand side the gain of population. Here 
Aij and Bij are the Einstein A and B coefficients for spontaneous and 
stimulated emission and absorption, J is the intensity of the radiation field 
averaged over all directions and integrated over the line profile, Cij is the 
collisional rate coefficient between levels i and j, and ricoi is the density of 
collision partners, primarily H2. The upward and downward rate coefficients 
are related through detailed balance by Cji = {gi/ gj)Cij exp{—hiy/kT). The 
values of Cij for collisions of different molecules with H2 J=0 have largely 
been derived from detailed quantum mechanical calculations (Green 1975, 
Flower 1990), and form one of the largest uncertainty in analyzing molecular 
excitation. Information on collisional rate coefficients at high temperatures 
>200 K is still scarce, and collisions with H2 J=l or higher are generally 
not taken into account explicitly. 

The critical density for a given transition is the density at which the 
rates for collisional processes become comparable to those for radiative 
processes, resulting in substantial population of level i. In formula, n*j,j^ = 
Aij/ J2j Cij- If the transition becomes optically thick, the critical density is 
lowered by 1/r, since part of the emitted photons are re-absorbed. Here r 
is the optical depth of the line. Note that higher frequency transitions have 
higher critical densities, since Aij ex z/^. Also, molecules with large dipole 
moments fi have high critical densities, since Aij oc fj?. Thus, by choosing 
the appropriate molecule and transition, a large range of densities can be 
probed. This is schematically indicated in Figure 7, which shows the range 
of physical parameters for which different line ratios are most sensitive. 

The contribution of the line photons to the excitation is taken into 
account by a variety of formalisms, including the large velocity gradient. 
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Figure 7. Examples of molecular line ratios which are used to constrain the density and 
temperature structure of YSO envelopes. 



Sobolev, microturbulent, or escape probability methods (e.g. 
Habing 1988, Osterbrock 1989 Appendix 2). 



Sobolev 1960, 



4.3. FROM OBSERVATIONS TO ABUNDANCES 

Figure 8 gives a schematic overview of the steps involved in the determi- 
nation of abundances from observations. Correcting for atmospheric effects 
and telescope losses yields the antenna temperature. This is the actual sky 
brightness averaged over the main forward beam of the telescope, or, for in- 
terferometric observations, the sky brightness averaged over the synthesized 
beam size after spatial filtering by the array baselines. 

Relating the antenna temperature to the underlying brightness distribu- 
tion on the sky requires assumptions about the source structure. Figure 8 
illustrates two extreme approaches to this problem. The 'homogeneous' 
method, shown on the left, uses a single set of physical conditions and a 
given size, or beam-filling factor, for the source, which transforms the beam- 
averaged antenna temperature into brightness temperature on the sky. The 
molecular excitation and opacity relate this brightness temperature to the 
column density of the molecule. The line opacity follows from the inten- 
sity ratio of lines of different isotopic varieties. The excitation can either 
be assumed to be in local thermodynamic equilibrium, valid if the density 
exceeds the critical density for the observed transition, or can be derived 
from statistical equilibrium and radiative transfer calculations as described 
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above once the physical parameters have been constrained. Repeating the 
procedure for an appropriate abundance standard hke CO yields the molec- 
ular abundance relative to H2. 

More complex source structures can be described with multiple 'homo- 
geneous' components, each with its own beam-filling factor. This approach 
works well if the number of different components is limited, for example, if 
the cloud is described as a collection of dense clumps embedded in a more 
diffuse interclump medium (e.g., Hogerheijde et al. 1995). When the beam 
is filled with a large number of different physical regimes, as is the case, 
e.g., for the power-law density and temperature distributions often used 
to describe YSO envelopes, the multi-component 'homogeneous' approach 
no longer suffices. The right-hand side of Figure 8 shows a more 'detailed' 
method. Starting with a fiducial source model, the molecular excitation 
and the radiative transfer throughout the cloud is calculated. Monte-Carlo 
techniques are often used in this step (e.g., Bernes 1979). The emergent sky 
brightness distribution is convolved with the telescope beam, or sampled at 
the interferometer's spatial frequencies. Matching the observed and model 
antenna temperature constrains the molecular abundance. 

The interpretation of interferometric observations, and the derivation 
of molecular abundances from them, requires that spatial filtering by the 
interferometer is taken into account explicitly. Even if the data sets of 
the molecule in question and that of a standard like C^^O contain the 
same spatial frequencies, both lines can have very different distributions 
of intensity over those frequencies if they trace different physical regimes. 
Different levels of flux are resolved out by the interferometer, resulting in 
abundance estimates deviating by factors of a few or more if not accounted 
for properly. 

5. Chemistry in pre-stellar cores 

In the following sections, an overview will be given of recent observations 
and models of star-forming regions at different evolutionary stages, using 
the material presented in §1-4 as general background information. The dis- 
cussion focuses on those species whose abundances are most affected at a 
particular evolutionary state. These molecules are often only minor com- 
ponents in terms of the overall composition of the gas or dust. Recent 
overviews of the major reservoirs of the principal elements C, N and O are 
given, e.g., in van Dishoeck et al. (1993) and van Dishoeck & Blake (1998). 
In brief, most of the oxygen is contained in solid silicates, oxides and ices, 
as well as gas-phase O, CO and, in warm regions, H2O. Most of the carbon 
is in some solid carbonaceous form, with the remainder in gas-phase CO. 
Nitrogen is mostly locked up in gas-phase N2 and N, but the amount in 
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Figure 8. The steps involved in deriving molecular abundances from spectral line data. 
The left-hand side shows the 'homogeneous' approach, the right-hand side the 'detailed' 
approach. 



solids is poorly constrained. 



5.1. TRANSLUCENT CLOUDS 



In order to study the effect of YSOs on the chemistry and test the basic 
chemical networks, the abundances in quiescent clouds prior to star for- 
mation need to be constrained. The best clouds for this purpose are the 
translucent and high-latitude clouds with visual extinctions of a few mag 
and densities of a few thousand cm~^. These clouds have been studied by 
optical absorption lines against reddened stars (e.g., Gredel et al. 1993), 
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by millimeter absorption lines against distant radio sources (e.g., Lucas &: 
Liszt 1997) and by millimeter emission observations (e.g., Turner 1996, 1998 
and references cited). No signs of star formation have been found in these 
regions. The chemistry is characterized by simple diatomic and triatomic 
molecules, radicals and ions. Detailed models by Turner et al. (1998) show 
that the observed abundances of most species are well reproduced by the 
ion-molecule networks. The main exceptions are formed by NH3, H2CO and 
H2S, whose formation is probably dominated by grain-surface chemistry. 

5.2. DARK CLOUD CORES 

More than 100 dark cores have been identified through optical extinction 
and molecular line studies, at least half of which have not yet formed any 
stars (Myers & Benson 1983, see chapter by Myers). Most chemical stud- 
ies have been directed toward TMC-1, a small elongated, dense ridge of 
~ 0.6 X 0.06 pc in Taurus with a mass of less than 40 M©. This clump 
shows a particularly rich chemistry, with a large chemical gradient across 
the core. NH3 and other 'late' molecules peak in the northern part, whereas 
long carbon-chain molecules such as HC7N and C4H are most abundant in 
the southern part (Olano et al. 1988, Pratap et al. 1997). High spectral 
and spatial resolution observations show that there are at least three differ- 
ent velocity components in this region on less than 10,000 AU scales, with 
different intensity ratios among the molecules (Langer et al. 1997). The 
abundances of the carbon-chains range from 10~^^ up to 10~^ and are sev- 
eral orders of magnitude larger than found in other dark clouds. Such large 
abundances are usually interpreted with pseudo time-dependent models at 
'early times' t ~ 10^ yr, assuming that the cloud evolved from a diffuse 
cloud phase in which carbon was initially in atomic form (e.g., Lee et al. 
1996a). 

In order to investigate whether TMC-1 is chemically peculiar, Suzuki 
et al. (1992) and Benson et al. (1998) performed a systematic study of 
a few characteristic molecules (C2S, HC3N, HC5N, N2H+ and NH3) in a 
large set of dark cores. The abundances of the carbon chains are found to 
correlate well with each other, but not with NH3 or N2H"'", consistent with 
the case of TMC-1. The observed C2S/NH3 abundance has been reproduced 
quantitatively in models which start from diffuse gas and form dense cores 
over a period of 10^ to 2 x 10^ yr (see Figure 9). In this scenario, TMC-IS 
is one of the youngest clouds. There are numerous other non-equilibrium 
processes, however, such as penetration of ultraviolet radiation, outflows 
or turbulence, which can lead to the breakdown of CO into C or C^ and 
therefore give the cloud core a chemically 'young' appearance. Thus, the 
'chemical age' does not necessarily measure the true age of the core, but 
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Figure 9. Observed CCS abundance versus the CCS/NHg abundance ratio in dark clouds 
with and without young stars. The sohd hue indicates the model abundance ratios for 
T—IQ K and n(H2)=10'' crn"^. The numbers along the line indicate the 'age' of the cloud 
in lO"" yr, assuming that the gas is initially atomic-carbon rich (from: Suzuki et al. 1992). 



indicates the time since some dynamical event produced an atomic-carbon 
rich phase. 

The model results are also sensitive to the [C]g/[0]g elemental abun- 
dance ratio in the gas. If this value is increased from the canonical ratio 
of 0.4 to higher values, improvement with observations for the complex 
molecules is obtained (Bergin et al. 1997, Terzieva &: Herbst 1998). Pratap 
et al. (1997) show that the observed chemical gradient in TMC-1 can be 
reproduced by changes in [C]g/[0]g along the ridge caused by a density 
gradient. 

This point is further illustrated by high spatial resolution observations of 
the L1498 pre-stellar core by Kuiper et al. (1996), which show a chemically 
differentiated structure with NH3 peaking in the inner region and C2S in the 
atomic-carbon rich (clumpy) outer part. Further chemical studies of such 
cores, especially those which have a more centrally concentrated density 
structure and appear to be on the verge of collapse (e.g., Ward-Thompson 
et al. 1994, Motte et al. 1998), are warranted. 

The ice mantles in quiescent dark clouds can be probed through the 
observation of field stars behind the clouds (e.g., Chiar et al. 1995, Whittet 
et al. 1998). H2O, CO and CO2 ices have been detected, which indicate 
that up to 40% of the heavy elements may be frozen out at densities of a 
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few xlO'^ cm~^ (Schutte 1999). Evidence for significant CO depletion deep 
inside dark cfouds {Ay > 20 mag) also comes from comparison of C^^O 
observations with extinction maps derived from infrared star counts (Lada 
et al. 1994, Kramer et al. 1999). 

5.3. IONIZATION FRACTION 

The ionization fraction x(e)=n(e)/n(H2) is an important parameter in the 
dynamical evolution of star-forming regions (see chapters by McKee, Shu, 
and Mouschovias) . In model calculations, x{e) drops from its high value of 
~ 10~^ at the edge of the core to ~ 10~^ in the center, and scales roughly 
with iC/nny^^ (Mihar 1990, Bergin k Langer 1997). At the edge, C+ is 
the main supplier of electrons, whereas deep inside metal ions such as Mg^ 
and Fe"^, as well as molecular ions such as Hg", HCO^ and HaO^ dominate. 
Recent observations of DCO^/HCO"'' in a sample of dark cores give 
ionization fractions of a few xlO^^, for an adopted cosmic ray ionization 
rate of 5 x 10^^^ s^^ (Plume et al. 1998, Williams et al. 1998). Surprisingly, 
no systematic trends between pre-stellar cores and cores with stars have 
been found, but this may be due to the large scale of the observations 
(-^10,000 AU). A somewhat lower ionization fraction of 10^^-10^^ is found 
for more massive cores (Bergin et al. 1999, de Boisanger et al. 1996). 

6. Chemistry in cold envelopes around YSOs 

6.1. MODELS 

In the initial stages of collapse, the density increases strongly whereas the 
temperature stays low, T ~ 10 K. The principal prediction of the mod- 
els is enhanced freeze-out of molecules onto the cold grains under these 
conditions. The only exceptions are H2, He, HJ^ and perhaps N2. Models 
appropriate for the cold outer envelopes have been made by Rawlings et 
al. (1992), Willacy et al. (1994), Bergin k Langer (1997) and Shalabiea k 
Greenberg (1995), using parametrized fits to the density profiles in simple 
collapse models such as that of Shu (1977) or Basu k Mouschovias (1994). 
They differ strongly in the adopted mechanisms that return molecules to 
the gas, ranging from none in Rawlings et al.'s to efficient desorption in Sha- 
labiea k Greenberg's. The results are also sensitive to the adopted binding 
energies on the icy mantles, especially whether the outer layer is H20-rich 
or CO-rich. The ions HCO^ and N2H'^ are predicted to be good tracers 
of cold envelopes, because their abundances remain high owing to the in- 
crease in the Hg" abundance when its main removal partners (CO, O, ...) 
are depleted. The use of HCO^ to trace the envelope structure has been 
demonstrated by Hogerheijde et al. (1997, 1998). 
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Figure 10. ISO-SWS spectrum of the deeply embedded massive YSO W 33A. Various 
absorption features due to silicate grain cores and icy mantles are indicated (from: Gibb 
et al. 1999). 



6.2. OBSERVATIONS 



Observationally, depletion in collapsing cores is very difficult to prove, be- 
cause every cloud has a 'skin' in which the abundances are close to normal 
(Mundy & McMullin 1997). Even if the skin amounts to only a few % of 
the total column density, its factor of 10-100 higher abundances can effec- 
tively mask any depletions deep inside. The dust obscuration in the deeply 
embedded class stage is still so high that the young stars are too weak at 
near- and mid-infrared wavelengths for direct observations of ices. Careful 
modeling of the line and dust emission appears the only way to probe the 
abundances. One of the best-studied cases is that of NGC 1333 IRAS4, 
where depletions of more than a factor of 10 have been inferred (Blake et 
al. 1995). This phase of high depletion appears short-lived, however, since 
for other class objects like Serpens SMM 1 much less freeze-out has been 
inferred (see §8.1). 

At the more evolved, but still deeply embedded class I stage, infrared 
absorption line studies of the ice mantles in the cold envelopes become 
possible, both for low- and high-mass objects (Whittet 1993, Tielens & 
Whittet 1997, Chiar et al. 1998, Teixeira et al. 1998). Ground-based and 
ISO studies allow a nearly complete inventory of these ices to be obtained 
for deeply embedded massive young stars such as NGC 7538 IRS9, W33 A 
and RAFGL 7009S (Whittet et al. 1996, d'Hendecourt et al. 1996, Ehren- 
freund et al. 1997) (see Figure 10). The molecules that have been identified 
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and their typical abundances are summarized in Table 3. Although only 
abundances down to ~ 0.5% of H20-ice (~ 10^^ with respect to H2) can 
be probed, the inferred mantle composition is remarkably simple, consist- 
ing of species resulting from the hydrogenation and oxidation of O, C, N 
and CO. The total column density of gas-phase CO is larger than that of 
solid CO integrated along the line of sight for all sources studied so far. 

7. Chemistry in w^arm envelopes around ne^vly- formed stars 

When a young star begins to heat the inner surrounding envelope by radia- 
tion and/or shocks above ~ 100 K, various distinct chemical regions occur, 
including the hot core region, the region of interaction of the outflow with 
the surrounding envelope, and the more extended radiatively-heated warm 
envelope (see Figure 11). 

7.1. WARM ENVELOPES 

Self-consistent models of the thermal balance, chemistry and radiative trans- 
fer in spherical envelopes heated by an internal source have been con- 
structed for low-mass YSO's by Ceccarelli et al. (1996) and for high-mass 
YSOs by Doty & Neufeld (1997). Calculations of the radial dust temper- 
ature distribution include those of Campbell et al. (1995), Kaufman et al. 
(1998), and van der Tak et al. (1999). For a ~ 10^ Lq source, the region 
where T^ust ~ ^gas > 90 K extends to ~ 10^^ cm, resulting in evaporation 
of H2O and CH3OH ices and an enhancement of the gas-phase abundances 
of these species by a factor of 100-1000 (see Figure 11). The chemistry in 
this inner hot core is discussed in §7.2. The region where T > 20 K ex- 
tends to ~ 10^^ cm. An increase in the temperature from 10 to 100 K does 
not have a significant effect on the gas phase chemistry, but does result in 
evaporation of volatile species like CO. 

Observational data on molecules in star-forming regions are scattered 
throughout the literature, using a variety of (sub-) millimeter telescopes 
with different beam sizes. Comprehensive single-dish studies of a number 
of molecules have been performed for only a few sources, including W 3 IRS5 
(Helmich & van Dishoeck 1997, see §8.2) and NGC 2264 IRSl (Schreyer 
et al. 1997). More detailed physical models are needed to disentangle the 
contributions of changing temperatures and densities in the envelope from 
possible radial abundance gradients. 

Important complementary information is obtained from high-resolution 
infrared observations. Mitchell et al. (1990) showed the presence of both 
cold (T < 60 K) and hot (T = 120-1000 K) gas along the lines of sight 
toward massive YSO's from ^^CO and ^^CO observations. The hot gas 
contains high abundances of C2H2, HCN and CH4 (Lacy et al. 1989, 1991; 
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Figure 11. Schematic illustration of the physical and chemical environment of massive 
YSOs. The variation in the structure of the ice mantles due to heating and thermal 
desorption is shown (based on Tielens et al. 1991, Williams 1993, van Dishoeck & Blake 
1998). 



Evans et al. 1991; Carr et al. 1995; Boogert et al. 1998; Lahuis &: van 
Dishoeck 1999), which are enhanced by at least 1-2 orders of magnitude 
compared with the colder envelope. Gas-phase CO2 has a surprisingly low 
abundance (van Dishoeck et al. 1996, Dartois et al. 1998, van Dishoeck 
1998b). Clear variations in the gas/solid ratios are seen for various objects, 
indicating the development of a hot core in the inner envelopes. 

7.2. HOT CORE REGIONS 

Hot cores are small (<0.1 pc), dense (n(H2)> 10^ cm~^) and warm (T sa 
200 K or higher) regions located close to massive young stars (see Walms- 
ley &: Schilke 1993, Ohishi 1997 for reviews). Observationally, hot cores are 
characterized by high abundances of fully hydrogenated molecules such as 
NH3, H2S, and H2O (e.g., Gensheimer et al. 1996), as well as saturated com- 
plex organic molecules hke CH3OH, CH3CN, CH3OCH3 and HCOOCH3 
(e.g., Hatchell et al. 1998). Examples of hot cores include the Orion hot core 
and compact ridge (e.g., Blake et al. 1987, Sutton et al. 1995), SgrB2(N) 
(e.g., Kuan & Snyder 1996), W 3(H20) (e.g., Helmich & van Dishoeck 1997) 
and objects near ultracompact H II regions such as G34. 3-1-0. 15 (Macdon- 
ald et al. 1996; see chapter by Churchwell). The observed abundances of 
the species can vary significantly from region to region. 

The rich chemistry in hot cores is thought to be driven by the evapo- 
ration of icy mantles due to heating near the young star. Charnley et al. 
(1992) showed that the observed abundances in the Orion hot cores can 
be reproduced if a mixture of simple ices containing H2O, CO, CH3OH, 
NH3 and/or HCN is evaporated into the warm gas. These molecules sub- 
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Figure 12. Models of the gas-phase chemical evolution during the 'hot core' phase for 
r=200 K. At t=0, molecules such as CH3OH, NH3 and H2S are evaporated from the 
grain mantles into the gas phase, where they drive a rapid chemistry leading to more 
complex organic molecules such as CH3OCH3, HCOOCH3 and CH3CN. The abundances 
peak after 10* — 10'' years (based on Charnley et al. 1992, Charnley 1997). 



sequently drive a rapid gas-phase chemistry leading to complex organic 
molecules (e.g., Charnley et al. 1992, Caselli et al. 1993, Charnley 1997). 
The abundances of the complex species peak after ~ 10^ yr, thus offering a 
chemical 'clock' since the time of formation of the hot core (see Figure 12). 
At later times, the normal ion-molecule chemistry takes over, resulting in 
destruction of the complex organic molecules after ~ 10^ yr. 

Since methanol is one of the most abundant ice mantle consitutents, it 
plays a crucial role in the hot core chemistry. Specifically, complex oxygen- 
containing organics can be formed by protonation of CH3OH by H^, HsO"*" 
or HCO^, followed by transfer of a -CH3 group to a neutral molecule, e.g., 
CH3OH+ + CH3OH -^ (CH3)20H+ + H2O. Dissociative recombination 
of the ion then produces di-methyl ether, CH3OCH3. The detailed chem- 
istry of interstellar alcohols leading to large esters and ethers has been 
discussed by Charnley et al. (1995). Nitrogen-containing organics can be 
formed through reactions of carbon ions with evaporated NH3, or through 
reactions of CN with evaporated C2H2. In general, the time scale for for- 
mation of nitrogen-rich molecules is somewhat longer than that of oxygen- 
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bearing species. 

The H2S/SO2 ratio may be a particularly sensitive clock (Charnley 1997, 
Hatchell et al. 1998). Most of the sulfur is thought to be in the form of H2S 
in the ices, although this has not yet been confirmed observationally. Upon 
evaporation, H2S is destroyed by reactions with atomic H to atomic S at 
temperatures of a few hundred K. S can subsequently react with OH and 
O2 to form SO and SO2 on a timescale of ~ 10^ yr. The observed high D/H 
ratios in hot core molecules (e.g., Millar & Hatchell 1998) largely reflect the 
high D/H ratios in the ices, which result from the cold dark cloud phase 
prior to star formation. 

7.3. OUTFLOW REGIONS 

The earliest stages of star formation are accompanied by powerful out- 
flows, which originate within the young star/accretion disk boundary re- 
gion. As the high velocity gas flows outward, it strikes the envelope and 
creates shocks, which alter the chemistry. Recent reviews of the physical 
and chemical structure of shocks can be found in Draine & McKee (1993), 
HoUenbach (1997) and Bachiller (1997, this volume). 

If the shocks have high velocities (>50 km s^^), the temperature in- 
creases to such high values (> 10^ K) that most molecules are collisionally 
dissociated. In addition, the ultraviolet radiation from such J-type shocks 
can dissociate molecules both within and ahead of the shock. The molecules 
reform slowly in a lengthy, warm zone in the wake of the shock (Neufeld &; 
Dalgarno 1989a,b). 

For lower shock speeds in clouds with low fractional ionization (xg < 
10~^), the shock is of C-type and the temperatures reach peak values of 
only 2000-3000 K, insufficient to dissociate the molecules. Reactions with 
energy barriers such as O -|- H2 and S -|- H2 rapidly drive the available 
gas-phase oxygen and sulfur into H2O and H2S. 

In addition to these gas-phase processes, grain cores and mantles are 
affected by shocks. The high-velocity J-shocks can cause destruction or 
thermal sputtering of the cores, resulting locally in much enhanced gas- 
phase abundances of Si, Fe, ... The lower velocity C-shocks can inject a 
variety of refractory and volatile species into the gas phase through non- 
thermal sputtering (e.g., Pineau des Forets & Flower 1997). The released Si 
can then react with OH and O2 to form SiO, one of the principal diagnostics 
of shocks (e.g., Schilke et al. 1997). 

Outflows have been studied observationally in a variety of molecules, 
including H2 near-infrared emission, CO and SiO. One of the best studied 
regions chemically is LI 157 (Bachiller & Perez Gutierrez 1997). Both re- 
fractory (SiO, SO, ...) and volatile mantle species (HON, H2CO, CH3OH, 
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TABLE 4. Physical regions of class YSOs traced by various observations 



Region Observation 



Bulk of the (cold) envelope Interferometer: Continuum on short spacings 

Single-dish: Lines tracing Tkin < 40 K 
(e.g., C^^O 2-1; C'^O 3-2; H(i3)C0+ 3-2, 4-3) 

Warm, inner regions of the envelope Interferometer: Continuum on long spacings 

Single-dish: Lines tracing ~ 100 K 
(e.g., i^CO 6-5; H2CO 322-221) 

Outflow: gradual entrainment (?) Interferometer: V or U shaped structures 

(e.g., HCO+, HCN) 
Single-dish: line wings (e.g., CO, HCN) 

Outflow: direct impact Single-dish and interferometer: Products of 

ice mantle evaporation and grain destruction 
(e.g., SiO, SO) 



...) are found to be enhanced in the gas by factors ranging from a few to 10^. 
Strong H2O lines have been observed in outflow regions from the ground 
(e.g., Cernicharo et al. 1994) and with the ISO satellite, indicating H2O 
abundances of a few times 10~^ (e.g., Liseau et al. 1996, Ceccarelli et al. 
1998, van Dishoeck et al. 1998). 

8. Examples 

8.1. AN INTERMEDIATE-MASS CLASS YSO IN SERPENS 

Recent work by Hogerheijde et al. (1999) of the class YSO Serpens SMM 1 
(FIRS 1) illustrates the use of molecular- line and continuum observations 
with single-dish and aperture-synthesis instruments to constrain the phys- 
ical and chemical conditions in the envelope of a low- to intermediate-mass 
protostar. The Serpens molecular cloud {d ^ 400 pc) is in the process of 
forming a loosely bound cluster of more than 50 stars (Eiroa & Casali 
1992). SMM 1 is one of the submillimeter-continuum class YSOs iden- 
tified by Casah et al. (1993) with Lboi = 77 L© and M^ = 0.7-3.9 Mq. 
Table 4 gives an overview of the different physical components traced by 
the various observations. 

The continuum observations are useful to constrain both the mass and 
the density structure of the envelope. The total mass is derived from the sin- 
gle dish data, whereas the density variation of the envelope is best obtained 
from the spatially resolved emission in the interferometer. The method is 
illustrated on the left-hand side of Figure 13. The dependence of the flux on 
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Figure 13. Left: Schematic illustration of single-dish versus interferometer observations. 
Right: Comparison of model visibility amplitudes at 2.7 mm to observations of SMM 1. 
The observed vector-averaged amplitudes are indicated by the filled symbols and their 
1(7 error bars. The dashed lines are models without a central point source, and density 
power-law slopes of —1.0 {lower curve in each panel), —2.0 [middle curve), and —3.0 
[upper curve). Solid lines are models with a point source flux of 0.13 Jy at 2.7 mm and 
a spectral slope of 2.0. 



projected baseline length for SMM 1 (right hand side of Figure 13) indicates 
a power-law density distribution; different prescriptions like, e.g., a Gaus- 
sian distribution, are excluded by the data. A power-law index of —2.0 it 0.5 
gives a good fit to the observations, which agrees well with models for pro- 
tostellar collapse (e.g., Shu 1977). The fluxes on long baselines suggest the 
presence of extra, unresolved emission, probably associated with the inner 
few hundred AU of the envelope where the temperature is likely to exceed 
the adopted T^ust of r""'^ distribution. The absolute temperature scale is 
~27 K at a characteristic radius of 1000 AU, constrained by the spectral 
energy distribution over millimeter to far-infrared wavelengths. The total 
envelope mass is 8.7 Mq, with an uncertainty of a factor of 2-3 due to the 
adopted dust emissivity (Ossenkopf & Henning 1994). The fact that the 
inferred envelope mass is larger than the stellar mass confirms SMM I's 
nature as a class YSO. 

Using this description for the physical structure of the envelope, the 
molecular abundances can be determined using models of the line emis- 
sion. Hogerheijde et al. (1999) adopt a Monte-Carlo technique to solve the 
radiative transfer and the molecular excitation throughout the envelope. 
In addition to the physical conditions, the molecular line calculations re- 
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Figure 14- Cleaned, naturally weighted images of molecular line emission observed with 
OVRO toward SMM 1. Contours are drawn at 3a intervals, of 1 K km s^^ for HCO^ 
and SO, 8 K km s^^ for H^^CO^, and 4 K km s~^ for SiO.The synthesized beam size is 
indicated in each panel. The dashed circle shows the primary beam size. The arrows in 
the HCO"'' panel indicate the position angle of the radio jet from Rodriguez et al. (1989) 
(from: Hogerheijde et al. 1999). 



quire knowledge of the velocity structure. Systematic velocity fields like 
infall and outflow can be included, but in the simplest analysis a turbu- 
lent line width of --^1.4 km s~^ independent of radius is adopted, based on 
optically thin C^^O 3-2 and H^^CO^ 3-2 lines. The neglect of systematic 
velocity fields influences the results for optically thick lines like ^^CO and 
HCO"''. Note that the lines themselves can in principle also constrain the 
density structure (e.g., Hogerheijde et al. 1997), but are less accurate for 
low- to intermediate-mass YSOs because of lack of spatial resolution in the 
single-dish data and of signal-to-noise in the interferometer observations. 

The envelope mass can also be constrained from the optically thin CO 
emission. If the 'standard' dark cloud CO/H2 abundance of 10^'^ is adopted, 
the model overestimates the emission in the C^^O 2-1 and C^^O 3-2 lines 
by a factor of 2. Apart from adjusting the CO abundance, a gas kinetic 
temperature of 0.6 times the dust temperature reconciles this discrepancy; 
^kin < ^dust IS expected because of imperfect thermal coupling between the 
dust and the gas, and line cooling of the gas. Alternatively, CO can freeze 
out onto dust grains when the temperature drops below its sublimation 
temperature of ~20 K (see Table 3), lowering the emission. A CO depletion 
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of a factor of 6 in the cold outer layers of the envelope provides a good fit to 
the C^^O 3-2 and C^^O 2-1 lines. These lines are particularly good tracers 
of depletion below 20 K, since the 2-1 lines, tracing T^in -^ 16 K gas, are 
stronger affected than the 3-2 lines, tracing Tkin ~ 30 K gas. The inferred 
depletion of a factor of ~6 for SMM 1 is smaller than found toward other 
class envelopes (e.g., 10-20 for NGC 1333 IRAS 4, Blake et al. 1995), 
suggesting that the evolutionary phase characterized by large depletion 
factors may be short-lived, or that local environment (e.g., heating of the 
envelope from the outside) influences the depletion. 

Most of the observed submillimeter transitions of other molecules trace 
material in excess of ~ 30 K, and the derived abundances do not depend 
critically on the exact value of Tkin or possible depletion below 20 K. Table 5 
lists the inferred abundances for the parameters of the envelope derived 
from the dust emission. Compared to other well-studied class YSOs such 
as IRAS 16293-2422 (van Dishoeck et al. 1995), the abundances of HCO+, 
HCN, and II2CO are very similar, while C3H2, CN, and IIC3N are enhanced 
by an order of magnitude in SMM 1. The values for SiO and SO are lower 
by a factor of 10, most likely because the 14"-19" single-dish beam of the 
SMM 1 observations does not include the tip of the outflow where these 
species peak (see Figure 14). 

Many but not all lines are well fit by this model. High-excitation lines 
like ^^CO 6-5 and II2CO 322-221 indicate the presence of an additional 
amount of warm, ~ 100 K, material. The estimated column density of this 
gas amounts to less than 1% of the envelope mass, but dominates in emis- 
sion of these highly excited lines. It is likely that this material corresponds 
to the inner few hundred AU of the envelope, where the temperature may 
exceed the adopted distribution. Some of the emission in the other lines 
may originate in this component as well. Table 5 includes the abundances 
derived under the assumption that all emission originates in this warm gas. 
These values are larger by typically an order of magnitude, illustrating that 
care needs to be taken to separate these two contributions. 

Figure 14 shows the OVRO images of the low-excitation lines of var- 
ious molecules. The emission in the line wings of optically thick lines of 
CO, HCO"*" and HCN is clearly associated with the bipolar outflow, as is 
that of SiO 2-1. HCO+ outlines the outflow cavity, possibly tracing the 
slow entrainment of material into the flow. SiO coincides with the axis of 
the outflow, likely revealing destruction of dust particles by the direct im- 
pact of the outflow on ambient material. To derive abundances from these 
interferometer observations, the modeling has to include spatial filtering. 
For example, only 0.25 Mq is traced by the optically thin C^^O 1-0 in- 
terferometer observations, but if the resolving-out of extended emission is 
taken into account, this value is consistent with the total envelope mass 
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TABLE 5. Derived molecular abundances in the envelope 
of Serpens SMM r 



Species 


Envelope 


Warm gas'' 


IRAS 16293-2442'= 


i^CO 


= l(-4) 


= l(-4) 




HCO+ 


l(-9) 


2(-8) 


2(-9) 


HON 


2(-9) 


5(-8) 


2(-9) 


H2CO 


8(-10) 


9(-9) 


7(-10) 


C3H2 


2(-10) 


3(-9) 


4(-ll) 


CN 


5(-9) 


3(-8) 


I(-IO) 


HC3N 


2(-10) 


9(-10) 


3(-ll) 


SiO 


l(-ll) 


I(-IO) 


I(-IO) 


so 


2(-10) 


2(-9) 


4(-9) 



" From Hogerheijde et al. (1999). 

'' Abundances derived under the assumption that all emission originates in warm 100 K 

gas. 

'^ From van Dishoeck et al. (1995). 



of 8.7 Mq. Firmer conclusions about abundance variations on small scales 
requires more realistic, >1D models of the physical structure, including a 
correct treatment of the velocity field. 

8.2. THE W 3 MASSIVE STAR-FORMING REGION 

The W 3 massive star-forming region at ~ 2.3 kpc provides an excellent 
opportunity to study the chemistry of massive YSOs (L ^ IO^Lq) at dif- 
ferent evolutionary stages originating from the same parent cloud. Helmich 
et al. (1994) and Helmich & van Dishoeck (1997) performed an unbiased 
single-dish submillimeter spectral survey of three YSOs at 15" (0.15 pc) 
resolution: IRS5, IRS4 and W 3(H2 0). Although such observations lack 
the spatial resolution of the data on nearby low-mass objects discussed 
above, they provide a useful global picture of the chemistry. In addition, 
the envelope mass associated with these objects is significantly larger than 
for low-mass YSOs, so that even minor species can be detected. 

Figure 15 summarizes the spectra toward the three objects in the 345 
GHz atmospheric window. Clear physical and chemical differences are found 
between the three sources, in spite of their similar luminosities. The beam- 
averaged densities are at least 10® cm^'^ and the temperatures range from 
~ 55 K for IRS4 to at least 220 K for W 3(H20). Toward W 3(IRS5), silicon- 
and sulfur-bearing molecules such as SiO and SO2 are prominent. This 
source has a powerful, massive outflow and is presumably the youngest of 
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Figure 15. Summary of the JCMT 335-365 GHz line survey of three massive 
star-forming regions in the W 3 molecular cloud. The spectra were constructed from 
the observed line intensities obtained in double-side-band mode. Strong lines in com- 
mon in the three spectra are labelled in the W 3 IRS4 spectrum only. Large physical 
and chemical differences are found between the three regions, which are attributed to 
different evolutionary stages (Based on Helmich & van Dishoeck 1997). 



the three sources, showing also strong ice absorption features in its infrared 
spectra. Toward W 3(H20), discovered by Turner & Welch (1984), organic 
molecules like CH3OH, CH3OCH3, and CH3OCHO are at least one order 
of magnitude more abundant, indicating that it is well into the hot core 
phase. Finally, only simple molecules and radicals are found toward IRS4. 
One possible interpretation is that this object is more evolved and has 
already broken free from the parent cloud, setting up an H II region and 
PDR on the back side of the cloud. 

More detailed physical models, which constrained by multi-line observa- 
tions of high-dipole molecules like CS, H2CO, HON and their isotopes, are 
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needed to infer reliable abundances. Recent analyses of the massive YSO 
GL 2591 suggest a power-law density distribution with an index of —1.0 to 
— 1.5, somewhat shallower than found for low-mass YSOs such as SMM 1 
(Carr et al. 1995, van der Tak et al. 1999). 

9. Chemistry in circumstellar disks 

Most of this chapter has focused on the chemistry in the envelopes of deeply- 
embedded YSOs. Young stars with ages 1 — 5 Myr have dispersed their 
envelopes but are still surrounded by circumstellar disks with masses of 
~ 10~^ — 10~^ Mq and sizes of --^100 AU, comparable to those inferred 
for our primitive solar nebula (see Beckwith & Sargent 1996 for a review). 
Imaging of the gas in disks around T Tauri and Herbig Ae stars has so 
far been limited to CO and its isotopes, for reasons of sensitivity. Obser- 
vations of other molecules are only just beginning. Single-dish surveys of 
various molecules have been carried out for a few objects, in particular 
DM Tau and GG Tau (Dutrey et al. 1997), TY Hya (Kastner et al. 1997) 
and LkCa 15, MWC 480 and HD 163296 (Thi et al., in preparation). Mil- 
limeter interferometers are now also capable of imaging the more abundant 
molecules in disks, showing interesting morphological differences (Qi et al., 
in preparation). 

The detection of species like HCO+, CN and HNC indicates that both 
ion-molecule, photon-dominated and gas-grain chemistry play a role. The 
amount of depletion is high because of the high densities (typically 10^ — 10^ 
cm~^ in the outer disk) and appears to be species-dependent, with the 
most volatile molecules remaining in the gas furthest from the central star. 
Chemical models which couple the gas-dust chemistry with the dynamical 
evolution as the material is transported inwards have been developed re- 
cently by Aikawa et al. (1997, 1998) and Willacy et al. (1998). At small 
disk radii, irradiation from the central star and stellar X-rays may also af- 
fect the chemistry. These studies are important first steps to determine the 
connection between interstellar and solar nebula processes in the formation 
of icy planetesimals such as comets and Kuiper-Belt objects. 

10. Concluding remarks 

Significant progress has been made in the last decade on models and obser- 
vations of the chemistry in star-forming regions. New (sub-)millimeter data 
at higher angular resolution and frequencies allow the dense and warm com- 
ponents of the YSO envelopes to be probed directly. Ground- and spaced- 
based infrared observations provide the first complete inventory of solid- 
state species along the line of sight and put constraints on gas/solid ratios 
of abundant species. The large variations in observed abundances illustrate 
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that the chemistry clearly responds to the enormous changes in temper- 
ature and density during star-formation. Molecules freeze out onto grains 
in the cold pre-stellar cores and outer envelopes surrounding YSOs, where 
grain-surface chemistry produces new species. They are returned to the gas 
by heating due to radiation from the young stars and by shocks due to 
outflows impacting on the envelope. A rich chemistry can ensue in the hot 
gas, leading to abundant complex organic molecules and driving much of 
the available oxygen into water. Thus, these molecules provide important 
chemical and temporal diagnostics of the YSO environment. 

The derivation of abundances from spectral line data is far from sim- 
ple, however, and invariably includes assumptions about the underlying 
source structure. This physical structure can be constrained by analyzing 
spatially resolved continuum data, as illustrated for the case of Serpens 
SMM 1, and by using line ratios of appropriate molecules to infer tempera- 
tures and densities. Different approaches, ranging from the "homogeneous" 
to the "detailed" method, offer optimal strategies for different cases. In- 
terferometric data require special care, because the instrument is a spatial 
filter which acts differently on lines tracing different physical conditions. 

The development of comprehensive gas-grain models has made it pos- 
sible to put the observational results into a coherent framework of the 
evolution of gas and solid species throughout the formation of protostars 
to their incorporation into circumstellar disks and eventually comets and 
planetesimals. Such models are rapidly progressing beyond the (pseudo) 
time-dependent method, and attempts are being made to couple the chem- 
istry with (multi-)dimensional hydrodynamics. Accurate information on the 
basic molecular processes entering such models remains a prerequisite. 

Future instrumentation, such as the next generation submillimeter ar- 
rays and air- or space-borne infrared telescopes such as SOFIA, SIRTF, 
FIRST and NGST, will allow these studies to be extended to much higher 
sensitivity and/or smaller scales of only a few AU. Together they should 
provide a much clearer view on the connection between matter found in in- 
terstellar clouds and in new planetary systems, and on the origin of organic 
matter found in primitive objects in our own solar system. 
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